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Testing the starburst/AGN connection with SWIRE 
X-ray/70/im sources 



0^ 
O 
O 
(N 

a 

(N 

o 

u 

6 



> 

oo 
m 

o 
a^ 
o 

> 



M. Trichas^*, A. Georgakakis^, M. Rowan-Robinson^, K. Nandra^, D. Clements^, 
M. Vaccari^ 

^Astrophysics Group, Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2BW, UK. 

^National Observatory of Athens, L Metaxa k V. Pavlou, Athens 15236, Greece 

'•^Department of Astronomy, University of Padova, Vicolo Osservatorio 5, 1-35122 Padua, Italy 



Accepted 23 June 2009 



ABSTRACT 

We explore the nature of X-ray sources with 70 /ini counterparts selected in the SWIRE 
fields ELAIS-Nl, Lockman Hole and Chandra Deep Field South, for which Chandra 
X-ray data are available. A total of 28 X-ray/70 /xm sources in the redshift inter- 
val 0.5 < z < 1.3 are selected. The X-ray luminosities and the shape of the X-ray 
spectra show that these sources are AGN. Modelling of the optical to far-infrared 
Spectral Energy Distribution indicates that most of them (27/28) have a strong star- 
burst component (> 5OM0yr~^) that dominates in the infrared. It is found that the 
X-ray and infrared luminosities of the sample sources are broadly correlated, consis- 
tent with a link between AGN activity and star-formation. Contrary to the predictions 
of some models for the co-evolution of AGN and galaxies, the X-ray/70 /im sources 
in the sample are not more obscured at X-ray wavelengths compared to the overall 
X-ray population. It is also found that the X-ray/70 /xm sources have lower specific 
star-formation rates compared to the general 70 /im population, consistent with AGN 
feedback moderating the star-formation in the host galaxies. 
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1 INTRODUCTION 

A major recent development in extragalactic astronomy is 
the discovery that a large fraction of the spheroids in the 
local Universe contain a super-massive black hole (Magor- 
rian et al. 1998). Moreover, the black hole mass is tightly 
correlated to the stellar mass of the host galaxy bulge (eg 
Ferrarese et al. 2000; Gebhardt et al. 2000), suggesting that 
the formation and evolution of galaxies and the build-up of 
the super-massive black holes at their centers are intercon- 
nected. The physical process behind this fundamental ob- 
servation is under debate, with AGN-driven feedback mech- 
anisms proposed as a plausible option. 

Recent simulations which include AGN feedback and postu- 
late that AGN form during galaxy mergers (Hopkins et al. 
2005; 2007) have received much attention recently because of 
their success in reproducing a range of observed properties of 
both AGN (e.g. duty-cycle, Nh distribution; Hopkins et al. 
2005) and gala^xies (e.g. bimodality of the colour-magnitude 
relation; Blanton 2006) . In this scenario galaxy interactions 
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trigger the AGN and also produce nuclear starbursts that 
both feed and obscure the central engine for most of its 
active lifetime. AGN-driven outflows (i.e. feedback) eventu- 
ally develop, which at later stages become strong enough 
to rapidly quench the star- formation, leaving behind a red, 
passive remnant and allowing the AGN to shine unobscured 
for a short period (~ 10* yr). 

A prediction of these simulations is that there should be an 
association between AGN and starburst events. Examples of 
such composite objects are rather isolated in the literature 
however (Cid Fernandes et al. 2001; Page et al. 2004; Alexan- 
der et al. 2005), particularly at z ~ 1, close to the peak 
of the AGN luminosity density in the Universe (Barger et 
al. 2005; Hasinger et al. 2005). This observational fact may 
point to a more complex AGN/starburst connection than 
that adopted in current models. For example, observations 
show that X-ray and/or optically-selected AGN are often 
found in red early- type bulge-dominated galaxies (Grogin 
et al. 2002; Pierce et al. 2006; Nandra et al. 2007) with evi- 
dence for star-formation events that have terminated in the 
recent past (Kaufi^mann et al. 2003 but see Silverman et al. 
2009). These observations are consistent with a time lag be- 
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Table 1. X-ray/70/^m sample 



object no. 


RA 


dec 




r b 




Nl 












1 


242.03664 


54.94373 


(0.862) 




12.70 


2 


242.14854 


54.39154 


1.300 






3 


242.30495 


53.90828 


0.9923 




11.87 


4 


242.45093 


54.60308 


0.8730 




12.01 


5 


242.55586 


55.11386 


(0.542) 




11.85 


6 


242.74052 


54.68533 


(1.270) 




12.45 


7 


242.80289 


55.13984 


1.280 




12.45 


8 


243.26450 


54.72134 


(1.138) 




12.34 


Lockman 












1 


161.13074 


59.09380 


1.312 




12.91 


2 


161.19244 


58.82482 


(1.280) 




12.29 


3 


161.01555 


58.98413 


(1.208) 




12.14 


4 


160.84180 


58.82229 


(1.051) 




12.01 


5 


160.97850 


58.98404 


(1.004) 




12.00 


6 


161.05058 


58.91622 


(0.845) 




11.94 


7 


161.96295 


58.86107 


(0.762) 




12.30 


8 


160.84120 


59.21511 


(0.754) 




11.79 


9 


160.85416 


59.28613 


(0.690) 




11.58 


10 


161.76088 


58.75598 


(0.644) 




11.70 


11 


160.72220 


59.17268 


(0.614) 


11.96 




12 


161.21414 


59.32958 


0.463 




11.33 


13 


161.03317 


58.74367 


0.555 




11.76 


14 


161.24121 


58.56060 


(0.905) 




12.11 


15 


161.27065 


58.99651 


0.264 




10.57 


CDFS 












1 


52.84087 


-27.85638 


(1.032) 




12.42 


2 


52.87523 


-27.93405 


(0.660) 




11.98 


3 


52.95019 


-27.80056 


(1.080) 




12.33 


4 


53.01069 


-28.05587 


(1.158) 




12.24 


5 


53.02176 


-28.07097 


(0.675) 




11.58 





T d 


r e 


type£ 


Av£. 




log Nl 


(A220) 


11 


.45 


11.20 


Sbc 




43.66 


23 


.56 




12 


.73 


12.82 


QSO 


0.15 


44.41 


22 


.46 




12 


.17 


12.32 


QSO 




45.06 


21 


.61 




10 


.41 


11.31 


Sbc 


0.2 


43.84 


23 


.67 


(A220) 


11 


.00 


11.10 


Sbc 




42.87 


22 


.89 




11 


.85 


11.45 


Sbc 




43.81 


21 


.71 




11 


.95 


11.75 


Sbc 




44.68 


23 


.60 




11 


.20 


11.61 


Sod 




43.58 


20 


.11 




12, 


.51 


12.61 


QSO 


0.10 


44.09 


22, 


.79 








11.35 


Sod 


0.3 


43.21 


19, 


.83 








11.14 


sb 


0.6 


43.18 


22 


.38 




10 


.44 


11.59 


Scd 




43.43 


19 


.83 








11.25 


Sod 


0.1 


43.07 


23 


.65 


(A220) 






11.64 


E 




43.43 


22 


.79 








11.60 


E 




42.81 


91 


.o / 




10, 


.50 


11.50 


Scd 




41.98 


21, 


.42 




11, 


.08 


10.98 


Scd 




43.35 


23, 


.43 


(A220) 






10.78 


Scd 


1.0 


43.35 


22, 


.38 








11.63 


Scd 


1.2 


42.94 


22, 


.37 




10, 


.93 


11.23 


sb 


0.2 


43.60 


22, 


.45 




11, 


.56 


11.10 


Sbc 


0.2 


42.25 


22, 


.54 


(A220) 


10, 


.44 


11.46 


sb 




42.99 


22, 


.56 








9.77 


Sab 




42.83 


23, 


.74 


(A220) 


12, 


.32 


11.02 


sb 




43.20 


23, 


.30 


(A220) 






11.53 


Scd 


0.2 


41.74 


19, 


.85 


(A220) 


11, 


.13 


10.33 


sb 




42.25 


19, 


.84 




11, 


.64 


11.06 


Scd 




43.60 


19, 


.89 




11, 


.32 


10.90 


Scd 




42.98 


22, 


.48 



Photometric redshifts are denoted in brackets 
^ cirrus luminosity 
^ starburst luminosity 

dust torus luminosity 
^ Bolometric optical luminosity 
/ Optical SED best fit 
9 extinction 

^ Hard band (2 — 10 kcV) X-ray Luminosity in erg s~^ 
* Intrinsic column density in cm^'^ 



tween the peak of the star-formation and the BH accretion 
in these systems. 

Alternatively, observational effects may bias our under- 
standing of the properties of AGN hosts. A major problem 
in assessing the level of star-formation in AGN is the dif- 
ficulty in decomposing the stellar from the AGN emission, 
especially in the case of dust enshrouded systems. Combin- 
ing X-ray data with longer far-IR wavelengths (e.g. 70/im) 
can provide a handle on this issue. The X-rays provide the 
most efficient way of identifying active BHs, least biased by 
obscuration, while common wisdom has it that the far-IR 
is a clean star-formation diagnostic. The sample, which is 
by no means complete, presented in this paper, has been 
compiled based on the above criteria; objects close to the 
peak of the AGN luminosity density which have both X- 
ray and 70^m detections, represents the best candidates for 
composite objects. Because of instrumental limitations (e.g. 
sensitivity, confusion), any sample selected at > 70/im, like 



the one proposed here, represents the brightest and rarest 
examples of an underlying population that future more sen- 
sitive surveys with Herschel will unveil. 



2 DATA 

The sample of X-ray sources with 70/im counterparts stud- 
ied in this paper is compiled from the SWIRE survey (Lons- 
dale et al. 2003), one of the largest Spitzer legacy pro- 
grammes. It covers a total area of 49deg^ spread in 6 dif- 
ferent fields. The observations have been carried out in all 
seven photometric bands available to Spitzer, the IRAC 3.6, 

4.5, 5.8 and 8 fim bands (Fazio et al. 2004) and the MIPS 
24, 70 and 160 /xm bands (Rieke et al. 2004). The typical 
5a sensitivities is about 3.7, 5.3, 48, 37.7fj,Jy in the IRAC 

3.6, 4.5, 5.8 and 8 /xm bands respectively. The correspond- 
ing 5(7 limits for the MIPS bands are 230fiJy at 24/im and 
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Figure 1. Photometric against spectroscopic resliift estimates 
for the X-ray/70/im sources. Photometric redshifts are estimated 
using the SED fitting method described in the text. Spectroscopic 
redshifts are from either our own follow-up program (Trichas et 
al. in prep) or from the literature. The diagonal line corresponds 

^phot ~ ^spec- 



20, 120mJy at 70 and 160 Atm (Surace et al., 2005). Optical 
observations in at least 3 photometric bands are available 
for most of the SWIRE area (>70%). The optical data com- 
bined with the SWIRE mid- and far-IR photometry have 
been used to fit the observed SEDs with a combination of 
galaxy and AGN templates (Babbedge et al. 2004; Rowan- 
Robinson et al. 2005) in order to study the nature of infrared 
selected populations and to estimate photometric redshifts 
for over 1 million sources (Rowan- Robinson et al. 2008) . De- 
tails about the latest SWIRE data release can be found in 
Surace et al. (2005). In this paper we have used multiwave- 
length data from three SWIRE fields that have X-ray data 
available from Chandra. These are the ELAIS-Nl, the Chan- 
dra Deep Field South (CDFS) and the Lockman field. The 
overlap between the Chandra X-ray observations and the 
SWIRE survey in these three fields is about 2.5 deg^. De- 
tails about individual fields are listed below. 



2.1 ELAIS-Nl 

The entire 9deg^ area of SWIRE ELAIS-Nl field has been 
observed at optical wavebands as part of the Wide Field 
Survey (WFS, McMahon et al., 2001) using the Wide Field 
Camera (WFC) on the Isaac Newton Telescope (INT). The 
survey consists of 600 s exposures in five bands: U, g', r', i' 
and Z to magnitude (Vega, 5a for a point-like object) lim- 
its of: 23.4, 24.9, 24, 23.2, 21.9 mag respectively. The overall 
photometric accuracy of the INT WFS survey is 2%. Fur- 
ther details are given in Babbedge et al. (2004) and Surace 
et al. (2005). 

X-ray coverage of ELAIS-Nl consists of two Chandra sur- 
veys. The first one is our own, and consists of 30x5ks expo- 
sures with the ACIS-I instrument. There is significant over- 
lap between adjacent ACIS-I pointings, thereby providing 
nearly even coverage over a circular area of about 1.5 deg^. 
This survey has not been fully presented in a previous publi- 
cation, except briefly in Georgakakis et al. (2008), and there- 
fore we will describe it in some detail here. The data were 
obtained in VFAINT mode during Cycle 6, between Decem- 




Figure 2. Optical to far-infrared SEDs in uf^, of all 5 X-ray 
sources with 70/im from CDFS. Solid curves show the total pre- 
dicted SED. Parameters for model fits given in Table 1. 

ber 2005 and March 2006. The data were analysed following 
the methodology of Laird et al. (2009). The CIAO v3.2 data 
analysis software was used for the data reduction. After cor- 
recting for hot pixels and cosmic ray afterglows using the 
AClS_RUN_HOTPlx task, the level 2 event file was produced 
using AClS_PROCESS_EVENTS. Flares were identified and re- 
moved using the automated method described by Laird et 
al. (2009). Individual event files although largely overlapping 
are treated separately to avoid problems arising from merg- 
ing regions with significantly different PSFs. Images are con- 
structed in four energy bands 0.5-7.0, 0.5-2.0, 2.0-7.0 and 
4.0-7.0 keV. Energy weighted exposure maps for each of the 
four energy bands were constructed using the merge_all 
task by adopting a power-law spectrum with spectral index 
r = 1.4 to estimate the energy weights. For each of the 
30 Chandra pointings of the survey, source catalogues were 
constructed in the energy bands above adopting a Poisson 
probability threshold of < 4 x 10~®, which are then merged 
into a master source list following the methodology of Laird 
et al. (2009). In the case of duplicate sources in the overlap 
regions of adjacent pointings we keep the detection with the 
smallest off-axis angle. The count rates in the 0.5-7.0, 0.5- 
2.0, 2.0-7.0 and 4.0-7.0 keV bands are converted to fluxes 
in the standard 0.5-10, 0.5-2, 2-10 and 5-10 keV bands re- 
spectively assuming a power-law X-ray spectrum with index 
r — 1.4 and Galactic absorption Nh = 2 x 10^" cm^^ appro- 
priate for the ELAIS-Nl fleld. The limiting fluxes are esti- 
mated 1.0x10"", 3.4x10"^^ and 1.5x10"" erg"i s"^ cm"^ 
in the 0.5-10, 0.5-2 and 2-10 keV bands respectively. A to- 
tal of 545 unique sources are detected over 1.5 deg^ in at 
least one of the four energy bands to the flux limits above. 
The second Chandra survey of the ELAIS-Nl is that pre- 
sented by Manners et al. (2003; 2004). This consists of a 
single 75 ks Chandra ACIS-I pointing, which lies in middle 
of our circular geometry shallow survey. Here we use the X- 
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ray source catalogue produced by these authors. They detect 
127 unique sources to flux levels of 2.3 x lO"^'^, 9.4 x lO"^*^ 
and 5.2 x 10"^^ erg-^ s"^ cm'^ in the 0.5-8, 0.5-2 and 2- 
8 keV bands respectively. 

2.2 Chandra Deep Field South 

The SWIRE CDFS was observed with the MOSAIC II cam- 
era on the 4m telescope at Cerro Tololo. Fifteen pointings 
covered ~ 4.5 deg^ in four filters, U, g, r, i to 5cr depths of 
24.5, 25.4, 25, 24 magnitudes (Vega) respectively with an 
additional 1.5 deg^ survey in the Z filter down to 23.3 mag- 
nitudes. There is also a deeper 0.33 deg'^ pointing in U, g, r, 
i to 5(7 depths of 25.2, 25.7, 25.5, 24.5 magnitudes respec- 
tively (Lonsdale et al., 2003). 

The X-ray data are from the Extended CDFS (ECDFS) and 
the IMs CDFS. We use the source catalogues of Lehmer et 
al. (2005) for the ECDFS and Giacconi et al. (2002) for the 
CDFS. The Extended Chandra Deep Field South consists of 
4 ACIS-I pointings covering 0.3 deg^ with an exposure time 
of 250 ks each. A total of 592 sources are detected by Lehmer 
et al. (2005) which reaches sensitivity limits of 1.1 x 10~^®, 
6.7 X 10"" ergs"^ cm"^ in the 0.5-2 and 2-8 keV bands re- 
spectively. The ECDFS is flanking the 1 Ms CDFS survey. 
Giacconi et al. (2002) have detected 316 X-ray sources down 
to flux limits of 5.5 x 10"", 4.5 x 10"" ergs"^ cm"^ in the 
0.5-2 and 2-10 keV bands respectively. 

2.3 Lockman Hole 

The Lockman Hole field was observed in the U, g', r', and i' 
bands with the MOSAIC Camera at the Kitt Peak National 
Observatory (KPNO) Mayall 4m Telescope, February 2002 
(g', r', and i') and January 2004 (U band). The scale of the 
Camera is 0.26" /pix and the field of view is 36' x 36'. The 
astrometric mapping of the optical MOSAIC data is good 
to less than 0.4" and the seeing varied between 0.9 and 1.4 
arcsec. Data reduction weis performed with the Cambridge 
Astronomical Survey Unit (CASU, Irwin & Lewis 2001) 
pipeline, following the procedures described in Babbedge et 
al (2004). Fluxes were measured within a 3" aperture (di- 
ameter) and corrected to total fluxes using growth curves. 
Typical 5 sigma magnitude limits are 24.1, 25.1, 24.4 and 
23.7 in U, g', r' and i' respectively (Vega), for point-like 
sources (Berta et al 2007a). 

The X-ray data are from the medium deep Chandra survey 
of Polletta et al (2006). It consists of 9 ACIS-I contiguous 
pointings of 70 ks each which cover an area of 0.6 deg^. The 
catalogue of Polletta et al (2006) has a total of 827 unique 
sources to the 0.3-8 keV flux limit of about 10"^^ erg s cm"^ . 
The total number of Chandra X-ray sources in the three 
fields is 2091. 



3 THE X-RAY/70 fxM SAMPLE 

In order to identify X-ray sources with 70 nm counterparts 
we have cross-correlated the SWIRE catalogue with the X- 
ray source lists described above following the same proce- 
dure discussed in Fadda et al. (2002), Franceschini et al. 
(2005) and Polletta et al. (2006). The total number of X- 
ray/SWIRE sources with 24 nm detections is 1368. In order 
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Figure 5. 2-10 keV X-ray luminosity as a function of total IR 
luminosity, Lju. Black crosses arc the X-ray/70 sources. The 
continuous line is the X-ray/IR luminosity relation for local star- 
forming galaxies from Ranalli et al. (2003). All X-ray/70 jum 
sources deviate from this Ranalli et al. relation in the sense that 
they are X-ray luminous for their Lm, indicating that AGN ac- 
tivity dominates the X-ray luminosity. 



to explore the possible link between AGN activity and star- 
formation at the epoch close to the peak of AGN luminosity 
density (e.g. Bargcr ct al. 2005), wo further restrict the sam- 
ple to the redshift interval 0.5 < 2; < 1.3 using either spectro- 
scopic (if available) or photometric redshift estimates. This 
resulted in a total of 28 X-ray /SWIRE sources with 70/jm 
detections. Given the surface density of X-ray and 70 nm 
sources within this redshift interval, the expected number of 
random associations is < 2. This sample consists of 5 CDFS 
sources, 1 from the 1 Ms survey and 4 from the ECDFS, 8 
are ELAIS-Nl sources, 7 from the 5ks survey and one from 
the EDXS and 15 sources are in Lockman. Table 1 gives the 
positions, redshifts (bracketed if a photometric redshift) and 
rest-frame luminosities for the different optical, infrared and 
X-ray components fitted to these sources. 

3.1 Multiwavelength Properties 

In order to classify the SEDs and estimate photometric red- 
shifts of the sources in our sample, the [/-band to 4.5 /xm 
photometric data are first fitted using a library of 14 tem- 
plates, 3 QSO, 1 Starburst and 10 galaxy templates (Rowan- 
Robinson et al. 2008). Photometric redshifts are estimated 
at this stop. At the longer Spitzer wavebands, the stellar 
contribution is first subtracted from the photometric data 
by extrapolating the best fitting galaxy template from the 
previous step. The corrected 3.6-170 nm data are then fitted 
with a mixture of four templates, cirrus (quiescent), M82 or 
Arp220 starbursts, and AGN dust torus, providing informa- 
tion on the nature of sources and the dominant mechanism 
responsible for the infrared emission, AGN vs stellar pro- 
cesses. The availability of 70 /um detections for all 28 sources 
and 160 /xm detections for 7 of them is an advantage, pro- 
viding additional leverage for elucidating the infrared prop- 
erties of the sample. 

For the sample of 28 X-ray/70/im sources the minimum 
number of optical/near-IR bands (up to 4.5 //m) used for 
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Figure 3. Optical to far-infraxed SEDs in vfi, of all 8 X-ray sources with 70/im from ELAIS-Nl. Solid curves show the total predicted 
SED. Parameters for model fits given in Table 1. 
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Figure 4. Optical to far-infrared SEDs in vfv of all 15 X-ray sources with 70/um from Lockman. Solid curves show the total predicted 
SED. Parameters for model fits given in Table 1. 



the photometric redshift calculation is 5. This is sufficient 
to estimate reliable photometric redshifts. This is demon- 
strated in Figure 1, which plots photometric against spec- 
troscopic redshift for a total of 12 sources in the sample with 
optical spectroscopy available cither from our own follow- 
up programs (Trichas et al. in prep.) or from the litera- 
ture. The agreement is good, with an estimated accuracy of 
5z/{l-\-Zspec) ~ 0.02 for the photometric redshift determina- 
tions of the X-ray/70/Lim sources. Broad-band SEDs for all 



28 sources arc given in Figures 2, 3, 4. The optical galaxy 
or AGN fits to the f/-band to 4.5 /xm data show that the 
photometric redshift solutions are plausible. All but two of 
the 28 sources require an M82 or A220 starburst component. 
The two exceptions arc Nl-2, which appears to have only an 
AGN dust torus component in the infrared, and Lockman- 
11, which seems to be better fitted with a quiescent ('cirrus') 
component. 19 of the 28 sources require an AGN dust torus 
component. Since we believe the X-ray emission of all but 
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Figure 6. IRAC/MIPS colour-colour plot. Solid line represents 
the predicted obscured AGN parameter space (Lacy et al., 2004). 
Red dots are SWIRE 70/im sources within the redshift interval 
0.5-1.3. Blue dots are the 28 X-ray/70/im sources. 



two of the 28 galaxies is due to an AGN (see section 4) , the 
remaining 7 galaxies must be cases with weak dust tori (see 
discussion by Rowan-Robinson et al. (2009) of the Lockman 
SWIRE-CLASX sample). 



4 RESULTS 

The X-ray and mid-IR properties of the X-ray/70 /xm 
sources indicate that most of them host an AGN. Figure 

5 shows that the X-ray /70/im sources are offset to higher 
X-ray luminosities compared to the Ranalli et al. (2003) 
Lx — Liu correlation for starbursts. Only two sources in 
the sample, Lockman-7 and CDFS-2, may possibly be X- 
ray starbursts, since they have Lx /Lir ratios only slightly 
higher than the Ranalli et al (2003) mean, and show no evi- 
dence for AGN activity (infrared dust torus, broad emission 
lines, QSO optical SED). Also, the obscuration corrected X- 
ray luminosities of 26 out of 28 sources in the sample are 
Lx{2 - lOkeV) > lO"^ erg s"^ cm~^ higher than what is 
typically attributed to star-formation (Georgakakis et al. 
2007). The X-ray spectra of 19 X-ray/70 /.tm sources are 
hard, consistent with column densities Nh ^ 10^^ cm~^ (as- 
suming a power-law intrinsic X-ray spectrum with F = 1.9), 
thereby indicating an obscured AGN. At infrared wave- 
lengths the X-ray/70 /im sources have luminosities Lm > 
10^^ L0 and therefore belong to the class of Luminous and 
Ultra-Luminous Infrared Galaxies (LIRGs, ULIRGs). The 
template fitting method of section 3 provides information 
on the emission mechanism(s) that are responsible for these 
bright luminosities. The mid-IR part of their SEDs include 
a significant contribution (up to 90 per cent at 8 /im) from 
a hot dust component associated with AGN activity. This is 
further demonstrated in Figure 6 which shows the position 
of the X-ray/70 fim sources on the mid-IR colour-colour plot 
of Lacy et al. (2004). The majority are in the region of the 
parameter space occupied by QSOs (Georgantopoulos et al. 
2007) and power-law selected AGN (Donley et al. 2008). 
However, in addition to the contribution from AGN heated 
dust, our SED modelling also suggests a star- formation com- 
ponent in the mid-IR in almost all cases. This component be- 




Figure 7. Bolomctric infrared luminosity, Lir, versus 0.1-3.0 fim 
optical luminosity. The small filled circles are SWIRE 24 fim 
QSOs (black cirrus, red M82, green A220, blue AGN dust torus, 
data from Rowan-Robinson et al 2008). The black crosses are the 
subset of these with 70 /im detections. The red crosses are X- 
ray/SWIRE sources with 24 fim detections (1368 objects). X-ray 
luminosities are calculated as in RR, Valytchanov and Nandra 
2009. The large black filled hexagrams are the subset of these 
with 70 ^m detections. 



comes increasingly important at longer infrared wavelengths 
and is essential to fit the the 70 and/or 160 nm emission of 
27 out of 28 sources in the sample. Quantitatively, accord- 
ing to our template fits, for 26 out of 28 sources in Table 
1, > 50 per cent of the Lir is from cool dust associated 
with star-formation. The estimated SFRs are > 50 Mq /yr, 
thereby indicating intense starburst events. In summary, the 
multiwavelength properties of the X-ray/70 pim sources in- 
dicate they are systems where both the stellar population 
and the SBH are growing at the same time. It is interesting 
to explore whether these two processes are related in the 
sample sources. 

For optically selected QSOs for example, Netzer et al. (2007) 
and Serjeant & Hatziminaoglou (2009) found that the 60/im 
luminosity, which estimates the current star-formation rate, 
scales with the 5100 A continuum luminosity, which is a 
proxy of the total AGN power, suggesting the two processes 
are linked. Figure 7 shows an equivalent plot, the bolometric 
infrared luminosity {Lm; 3 — 1000 ^m) to approximate the 
luminosity of the current star-formation rate versus the op- 
tical bolometric luminosity, Lopt (0.1 — 3 /^m) for all SWIRE 
QSOs with 24 fim detections (data from Rowan- Robinson 
et al 2008), as a proxy of the AGN power (see figure caption 
for details on the estimation of luminosities). X-ray/70 /xm 
sources and SWIRE QSOs show the same broad correlation 
between Lm and Lopt suggesting an association between the 
formation of stars in these galaxies and the growth of the 
SBH at their centres. We note however, that the correla- 
tion in Figure 7 has substantial scatter, i.e. at a given Lopt 
the LiH spans about 2dex. We have also shown the corre- 
sponding distribution for the whole sample of SWIRE/X-ray 
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Figure 8. Specific star formation rate in Gyr'^ versus stellar 
mass in Mq yr~^ for SWIRE galaxies, colour-coded by opti- 
cal template type, the TO/^m/X-ray sources are overplotted as 
large black filled hexagrams. Apart from two objects (CDFS-1 
and CDFS-3) at the top left, there is a tendency for our sample 
to have somewhat high stellar masses (> 2xlO^'^M0) and some- 
what lower specific star formation rates {<3Gyr~^ ). Both the 
two low stellar mass objects in our sample have dust tori. 



sources with 24 fim associations (red crosses in Fig 7), with 
X-ray luminosities first converted to a bolometric luminos- 
ity assuming a bolometric correction to the hard X-ray lu- 
minosity of 27, and then to Lopt by dividing the bolometric 
luminosity by 2 (Rowan- Robinson et al 2008). X-ray galax- 
ies with very high values of Lab /Lopt are likely to be X-ray 
starbursts. 

In contrast with the 70 /im QSOs, optically selected samples 
like the PG QSO samples of Netzer et al (2007) and Lutz 
et al (2008) are biased towards lower values of LsbjLopt- 
Optical selection favours objects with higher Lopt, 70 ptm 
selection favours higher Lir- The 70/xm/X-ray sample (filled 
hexagrams in Fig 7) seem to cover the whole range of Lapt 
and Lir 

Figure 8 shows the specific star-formation rate versus stel- 
lar mass for the entire SWIRE 70 \ixa population, in the 
redshift interval 0.5-1.5, colour-coded by optical template 
type, overplotted by our sample. Apart from two objects, 
CDFS-1 and CDFS-3 with low stellar masses and high spe- 
cific star-formation rates, the 70/im/X-ray population is at 
the lower end of the specific star-formation rate distribution 
for 0.5 < z < 1.5, and tends to have a higher stellar masses 
than the general population, i.e. they avoid the most intense 
starbursts of similar stellar mass. 

It is also interesting that a large fraction of the X-ray /70 ^m 
sample (19/28) are obscured by column densities Nh > 
10^^ cm~^. Figure 9 compares the HR distribution of the X- 
ray/70/im sources with that of the overall X-ray population 
with SWIRE photometric redshifts in the interval 0.5-1.5. A 
Kolmogorov-Smirnov tests shows that the likelihood of the 
observed differences if the two samples are drawn from the 
same parent population is about 3 per cent. Therefore there 
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Figure 9. Hardness ratio distribution for X-ray/70 /xm sources 
and the entire X-ray source population identified with SWIRE 
sources in the redshift interval. The hardness ratio is defined as 
HR=(H-S)/(H+S) where H is the hard band (2-7keV) count rate 
and S is the soft band (0.5-2 kev) count rate. The KS test shows 
that the probability that the distribution are drwan from the same 
parent samples is 3%, i.e. there are no statisitically significant 
(<2(t) differences between the two populations. 



is no strong evidence (3 per cent or about 2a in the case 
of a normal distribution) that the X-ray/70 /im sources are 
more obscured than the overall X-ray population. 



5 DISCUSSION 

Chandra X-ray data within the SWIRE survey footprint are 
used to select a total of 28 X-ray sources with detections at 
70 /im, with additional 160/im for 25% of them, in the red- 
shift interval 0.5-1.3 to study the properties of galaxies that 
are best candidates for systems where the formation of the 
stellar population is taking place at the same time as the 
growth of the SBH at their centres. Study of the multiwave- 
length f/-band to far-IR SED of these systems indeed shows 
that the majority of the sample harbors powerful starbursts 
(27/28; > 50Modot/yr), while the X-ray data indicate AGN 
activity in at least 26/28 sources. 

The X-ray spectra of 19/28 sources are consistent with mod- 
erate obscuration, Nh > lO'^^cm"^. However, the obscured 
AGN fraction in the sample is similar to that of the overall 
X-ray AGN population in the same redshift interval. This 
contradicts some models for the co-evolution of AGN and 
galaxies in which the central engine is obscured by dust and 
gas clouds associated with circumnuclear starbursts (e.g. 
Fabian et al. 1999; Hopkins et al. 2006; Ballantyne 2008). Al- 
ternatively, our finding may indicate that a substantial frac- 
tion of the AGN has ongoing star- formation and that the X- 
ray/70 /im might not be very different in their star-formation 
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properties from the bulk of the X-ray AGN population. Sil- 
verman et al. (2009) for example, used the [Oil] 3727 line 
as proxy of the star-formation rate in X-ray selected AGN 
at a fs 1 and showed that the majority show evidence for 

star-formation. 

We also find that X-ray/70 /xm avoid the most intense star- 
bursts at a given stellar mass. This is consistent with models 
where outflows from the central engine moderate the star- 
formation by sweeping away the cold gas reservoirs of the 
host galaxy (e.g. Hopkins et al. 2006). This result also con- 
tradicts Silverman et al. (2009) who finds that X-ray se- 
lected AGN in the COSMOS field have, on average, higher 
star-formation rate compared to non X-ray sources in that 
field of similar stellar mass. 



6 CONCLUSION 

SWIRE infrared and optical data have been combined with 
Chandra X-ray observations in the ELAIS-Nl, Lockman 
Hole and CDF-S to study the nature of the X-ray/70 /xm 
sources in the redshift interval 0.5-1.3. The majority of 
these systems are starburst/AGN composites thereby pro- 
viding the opportunity to explore the relation between SBH 
growth and galaxy formation. The X-ray/70 fiui sources fol- 
low the broad correlation between star-formation rate (ap- 
proximated by the L/h) and SBH accretion rate (approxi- 
mated by Lopt) for broad line QSOs, suggesting that these 
two processes are linked. We do not find evidence that the 
X-ray/70 /im arc more obscured than the overall X-ray pop- 
ulation, contrary to the predictions of some AGN/starburst 
co-evolution models. Also, the X-ray/70 /itm population ap- 
pears to avoid the most intense starbusts for their stellar 
mass. This may indicate AGN driven moderation of the star- 
formation in the host galaxy. 
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